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We explore the problems in calculating the decay constant for heavy-light mesons using the quenched Wilson 
quark action for both heavy and light quarks. We find that the continuum limit of the decay constant is reasonably 
converged among different prescriptions after the continuum limit is taken. A number of technical problems 
associated with prescriptions are also addressed. 



1. Introduction 

Extrapolation to the continuum limit is an in- 
dispensable procedure to obtain physics from lat- 
tice calculations. For the system including heavy 
quarks (Q) this extrapolation might no longer be 
straightforward, since 0(mQa) errors can uncon- 
trollably be large. In this work we numerically 
explore the problem of the continuum limit for 
the decay constant of heavy-light mesons using 
the Wilson heavy quark action in the quenched 
approximation. 

2. Simulation 

Tabic |l| summarizes parameters we used in our 
simulations carried out on VPP-500/80 at KEK. 
The improvement we have made since Lattice 
95 [0 is the adoption of a smeared source for heavy 
quarks to suppress contaminations from excited 
states, and an addition of a run at j3 = 5.9. We 
use as a smearing function the wave function of 
heavy-light mesons as measured in the simulation. 
The size of lattice is chosen so that physical vol- 
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Table 1 

Simulation Parameters 
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"size I6 ]i x 40 24 3 x 64 32 3 x 80 

#conf. 150 100 100 

q^_ 1P (GeV) 1.84(03) 2.54(07) 3.36(11) 

ume is roughly constant against j3. We take the 
physical scale of the lattice estimated from the 
charmonium spin-averaged 1S-1P mass splitting. 

3. Kinetic versus pole mass 

The effective Hamiltonian for Wilson quark 
with a large mass rriQa takes the formal, 
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where all mass parameters may differ by 0(mQa). 

In Figure [l]we plot mki n a—m po i e a as a function 
of the bare quark mass for heavy-light and heavy- 
heavy mesons, where rriki n a is extracted from 
measured dispersion relations. We expect the two 
results to agree, since the binding energy should 
cancel in the mass difference. The data, however, 
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Figure 1. rriki n a — m po i e a for heavy-light (filled 
symbols) and (rriki n a — m po i e a)/2 for heavy-heavy 
(open symbols) mesons. Dashed curve represents 
tree level value for free Wilson quark. 

show that they do not, implying that we can- 
not determine the kinetic mass consistently from 
heavy- light and heavy-heavy systems 0. Even 
worse, as /? increases, m,ki n a—m vo i e a undershoots 
and becomes negative in a region rriQa<0.7. 

We plot in Figure || 2m,Qq — rngQ , with mqq 
mass of a heavy-light meson and tuqq for a 
heavy-heavy meson. This mass difference as ob- 
tained with the two mass definitions is com- 
pared with experimental values (asterisks). The 
pole mass gives values close to the experiment, 
while the results with the kinetic mass definition 
stays far away from the experiment. Whereas we 
naively expect that rrikin is the most adequate 
mass parameter of the system, as the kinetic en- 
ergy term dominantly controls the dynamics of 
the system, our results show that the kinetic mass 
measured from simulations suffers from the effects 
arising from higher order terms of 1/mq which 
receive large finite lattice corrections [Q. For this 
reason we cannot conclude that the use of the 
kinetic mass is superior to the pole mass. 

4. Heavy-Light Decay Constant 

We study the problem of the continuum limit 
by taking a number of alternative prescriptions 
for the meson mass parameter and the wave func- 
tion normalization. (IA) The simplest choice is 
use of the pole mass and the naive normaliza- 
tion (2i^Q) 1 / 2 . (IB) Alternatively, we may take 
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Figure 2. 2m,Qq — mqQ calculated with pole mass 
(filled symbols) and kinetic mass (open symbols) 
as compared with experiment(asterisks). 

(8K C - QKq) 1 / 2 = {2K Q ) 1 / 2 e m ^ rk)a ^ as the 
wave function normalization (referred to as the 
KLM normalization). (HA) The next choice is 
use of the kinetic mass obtained from simula- 
tions together with the KLM normalization. In 
this prescription O(toqo) corrections are removed 
from leading terms of 1/rriQ. (IIB) In order to 
avoid the pathology we have seen with the ki- 
netic mass, we may consider yet another defini- 
tion of the 'kinetic mass' estimated from the pole 
mass using the tree level relation^), m kiv 
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while keep- 



(meson) . / (quark) (quark) \ 

""pole + V m kin ~ m pole 

ing the KLM normalization. The tadpole im- 
proved perturbative Z-factor is used throughout 
our analysis. 

Figure |^ shows the scale invariant quantity 
$p(to p ) = [a s (mp)/a s (m B )] 2/ P° fp^/mp as a 
function of 1/mp for cases (IA) and (IIB). The 
naive normalization gives significantly smaller 
values than the KLM normalization does for 
heavy quarks; as is well known fpyjmp even de- 
creases towards a large quark mass. This con- 
trasts to case (IIB), in which it increases with the 
meson mass and is smoothly extrapolated to the 
value estimated in the static limit The method 
(IB) also approaches the static limit, the differ- 
ence from (IIB) being an 0(e~ mp ) finite mass 
correction rather than 0(l/mp)Q. 

Extrapolation to the continuum limit is il- 
lustrated in Figure ^ for fp> 3 assuming linear 
0(rriQa) behavior. Roughly speaking, all pre- 
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Figure 3. $p(top) versus 1/mp with the pre- 
scription IA (open symbols) and with IIB (filled 
symbols). Results with the static approximation 
are also shown at l/rap=0. 

scriptions give a reasonably convergent answer, 
and the variation from one to another is an 
0((mQa) 2 ) effect. At a more precise look, how- 
ever, the result for (IIA) is sizably deviated from 
the others. This is apparently caused by a patho- 
logical behavior of mki n a — m po i e a in the charm 
mass region as discussed above. We note that the 
simplest prescription (IA) gives a result in agree- 
ment with the others in spite of a large extrapo- 
lation. The difference from (IB) by an exponen- 

tial factor e m p° le a/ in the KLM normalization 
does not invalidate the linear extrapolation for 
this case since 0{(mQa) 2 ) terms are small for c- 
quark. 

The results for fg a are also well convergent in 
the continuum limit (Figure ||) . Although the re- 
maining systematic uncertainty cannot be prop- 
erly estimated for the prescription using the pole 
mass for rriQa is larger than unity for the 6-quark, 
we expect that the remaining errors are reason- 
ably small for (IIB). 

After continuum extrapolation we obtain 
f D = 202(8)tii MeV f Ds = 216(6)±g MeV 
fs = 179(11)1?! MeV f Bs = 197(7)t° 5 MeV 
where we take the value with prescription (IIB) 
for our central values. The first error is statistical 
and the second represents a spread over four pre- 
scriptions, which gives an estimate of 0((toqo,) 2 ). 
Finally we emphasize that convergence of differ- 
ent extrapolations to the continuum limit both 
from above and below bolsters the reliability of 
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Figure 4. Continuum extrapolation of fn s with 
four methods: I A (open circles), IB (open dia- 
monds), IIA (crosses) and IIB (filled circles) 
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Figure 5. Same as Figure |j for fs,- 

our results, at least for charmed mesons. We 
expect that the errors for bottom mesons are at 
most within the quoted errors. 
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